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ABSTRACT 

We determine the spatial distribution function of galaxies from a wide range of samples in the 2MASS 
survey. The results agree very well with the form of the distribution predicted by the theory of cos¬ 
mological gravitational many-body galaxy clustering. On large scales we find a value of the clustering 
parameter b — 0.867±0.026, in agreement with b — 0.83 ±0.05 found previously for the Pisces-Perseus 
supercluster. We measure b{9) as a function of scale, since this is a powerful test of the applicability of 
computer simulations. The results suggest that when galaxies clustered they were usually surrounded 
by individual, rather than by communal haloes. 

Subject headings: large-scale structure of universe — dark matter — gravitation — infrared: galaxies 
— galaxies: clusters: general — galaxies: statistics 


1. INTRODUCTION 

The spatial locations of galaxies may be described 
by many different statistics. These include distribution 
functions, low order correlation functions, multi-fractal 
dimensions, topological genus, power spectra, spherical 
harmonics, multipoles, minimal spanning trees, perco¬ 
lation and moments. All these statistics are related 
and each emphasizes different aspects of the distribu¬ 
tion ('e.g.. lSaslawll200(ll . Here we determine the observed 
galaxy distribution function from the 2MASS Catalogs 
with considerably greater precision than has previously 
been possible from smaller catalogs, and compare it with 
theoretical expectations. 

The galaxy distribution function contains informa¬ 
tion about correlations to all orders, about voids and 
underdense regions, and about clustering over a very 
wide range of intensities, sizes and shapes. Moreover 
it can be calculated analytically for cosmological grav¬ 
itational many -body cluster ing in both linear and non¬ 
linear regimes i|Saslawll200n() . 

These theoretical calculations assume, consistent with 
direct N-body simulations, that galaxy clustering evolves 
through a sequence of quasi-equilibrium states. This oc¬ 
curs because in any overdense region the local dynamical 
timescale is shorter than the global timescale for evolu¬ 
tion of average quantities, such as the overall density, 
velocity dispersion, and the ratio of gravitational poten¬ 
tial correlation energy to the kinetic energy of peculiar 
velocities. The global timescale for these average ensem¬ 
ble properties to change, becomes even greater than the 
Hubble time, R{t)/R, as local partially virialized struc¬ 
tures form. The Hubble expansion exactly cancels the 
long-range gravitational many-body mean held, not only 
in the usual Einstein-Friedmann models, but also in those 
models with a cosmological constant and quintessence. 
For the cosmological constant and quintessence models, 
one can straightforwardly extend the earlier derivations 

^ Operated by Associated Universities, Inc., under cooperative 
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llSaslawll200nl: iSaslaw & Fanglir996D that applied to the 
Einstein-Friedmann case. 

Under these conditions, the clustering of galaxies, 
each of which may be surrounded by its own dark 
matter halo, occurs mainly through their mutual grav¬ 
itational interactions. This differs from those cold 
dark matter models that are dominated by very large 
haloes containing galaxies moving mainly as test par¬ 
ticles within the gravitat ional held of these haloe s 
(e.g, the massive haloes in iFukiishige fc Makinoll2?)0^ . 
When the mutual gravitational interactions of indi¬ 
vidual galaxies dominate, clustering can be described 
by quasi-equilibrium thermodynamics llSaslawl 120001: 
ISaslaw fc FantdllQOa fSaslaw fc Hamiltonlll984H and sta¬ 

tistic al mechanics ll Ahmad et al.ll2?)?)2: iLeong fc Saslawl 
l2?iol . This theory applies for a wide range of initial con¬ 
ditions, including i nitial powe r law perturbation spectra 
with — 1 ^ n < 1 ijltohlll 9^ . Characterising the com¬ 
plete range of initial conditions that form a basin of at¬ 
traction for this theory remains an important unsolved 
problem. 

Anal ysis of the cosmological ma ny-body problem pre¬ 
dicted llSaslaw fc HamiltonI Il984^ that its distribution 
function fv{N), i.e. the probability that a randomly 
placed volume V in space contains N galaxies, has the 
form 


fyi^N) = -b)+ 

( 1 ) 

where 

N = riV (2) 


is the expected number for an average number density n. 
The strength of clustering is measured by 


^ W 2TTGmSn 
~ ~ 3T 


/ Gn,T,r)-r‘^dr (3) 
<v r 


for a spherical volume where W is the gravitational cor¬ 
relation energy, K is the kinetic energy of galaxy peculiar 
velocities indicating an ensemble average temperature T, 
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the average galaxy mass is m, and ^ is the two-galaxy 
correlation function. Since the volume V may be a cone 
with apex at the observer which projects the galaxies in 
its volume into a cell of area A on the sky, the same 
form of the distribution function in equation o applies 
to counts of galaxies in two-dimensional cells on the sky. 
This form of equation © assumes that evolution is negli¬ 
gible within the volume, which is reasonable for a narrow 
redshift band as in the 2MASS catalog (0 < z < 0.1). 
For volumes covering wide re dshift bands, ev o lution may 
be included as describ ed in iFang fc Saslawl lll997ll and 
iSaslaw fc Edead l)2000D . 

Although the form of equation o is independent of 
the shape and size of randomly located volumes which 
are not pathological (e.g., not fractal or chosen to avoid 
galaxies), the value of b will depend on volume and shape. 
This was found in early N-body simulations llltoh et al.l 
and explicitl y calculated for square p rojected cones 
using equation © C ./ahav fc Saslawll992ll . Comparisons 
of different values of N and b for cells of the same shape 
and size, but for different subsamples of galaxies, can be 
used to explore possible biases in these subsamples. If 
the bias is known a priori, it can be incorporated into 
equations JTL and © as described i n detail previously 
llLahav fc ^slawlllOT2t ISaslawl 1200011 . In principle, an 
unknown bias contained in subsamples can also be recov¬ 
ered from the results, but its uniqueness is more difficult 
to establish. _ 

The value of N can be determined directly from the 
catalog and the value of b follows for equation o from 
the observed variance of counts in cells of any particular 
volume and shape: 


{(ml) = 


N 


{i-b{v)y 


( 4 ) 


Therefore in principle the theory contains no free param¬ 
eters for comparing equation o with the observations. 
We can also fit N and b in equation CJ directly to the 
observed distribution function, and see if the slight dif¬ 
ference between these two methods of determining N and 
b provides useful information. 

After deriving equation o, its predicted form was 
found to agree well with counts in cells and void prob¬ 
abilities in several galaxy surv eys including the Zwicky 
Catalo g llSaslaw fc Crandl9^ . the UGC and ESO Cat- 
alogs llLa hay fc Saslawl I1992D and the IRAS Catalog 
l|Sheth er^bllTOT4j) ■ al l projected onto th e sky, as well as 
for the SSRS Catalog ( Fan£_^,Zou|199j)^ and the Pi^ es- 
Perseus Supercluster ( Saslaw fc Haaue-CoDilahlll998D in 
three dimensions. All previous comparisons have in¬ 
volved catalogs containing at most several tens of thou¬ 
sands of galaxies. The new catalogs becoming available 
have about 10-100 times as many galaxies, with differ¬ 
ent levels of sky coverage, homogeneity, and information 
about galaxy properties. 

These new larger catalogs are important because their 
improved statistics may reveal departures from equa¬ 
tion o resulting from different distributions of galax¬ 
ies and dark matter, merging, mass segregation, or other 
environmental effects. Section |21 describes the 2MASS 
Catalog and our method of analysis, Section © gives our 
results, and Section © describes some of their implica¬ 
tions. 


2. THE 2MASS CATALOG AND ITS ANALYSIS 

The Two Micron All Sky Survey employed telescopes 
in the Northern and Southern hemispheres to observe in 
the three infrared bands J (1.11 — 1.36/rm), H (1.50 — 
l.SO/im) and Ars(2.00 — 2.32 /j,to). Most of the sources 
are points, but 1,647,599 are resolved with respect to the 
observed 2MASS point spread function and they consti¬ 
tute the Extended Source Catalog (XSC). About 97% of 
these extended sources are galaxies^. Examination of the 
Abell 262 cluster suggests that most of the 2MASS galax¬ 
ies are local with z < 0.1 (Skrutskie, private communi¬ 
cation). In addition to the local population, the Catalog 
contains more distant very luminous galaxies and about 
0.4% of the catalog may be quasars with redshifts < 5 
(IBarkhouse fc 200111 . The low redshifts for nearly 
all the sources indicate that evolutionary corrections will 
not be significant for our statistics using this sample. 

To sele ct our sampl e s, we adopted Kron magnitudes 
following iMaller et al.1 ll2005ll and dereddened t hem us¬ 
ing the reddening maos of iSchlegel et alJ (119981) assum¬ 
ing A(J) = 0.82E{B - V), A{H) = QA8EfB - V), and 
A{K,) = 0.28E{B - E) dMa iewski et al ]| 2003 ) 3 . We 
required all galaxies brighter t han Ar.„ n = 12 to have 
Jo — Ksfl between 0.7 and 1.4 llMaller et alJl200,^ and 
all candidate galaxies to have star-galaxy separation in¬ 
dices typical of galaxies (Cscore < 1-4 and gscore < 1-4) to 
avoid confusion with nearby galactic sources. We also re¬ 
quire that a candidate is not listecPas a possible anomaly 
(e.g. artifact, cluster/galaxy piece. Galactic object or 
star) and that if a candidate is among the 25% visu¬ 
ally verified sources, it is not indicated as non-extended 
{vc 7 ^ 2). To minimize the relative contribution of stars 
in the samples, we usually remove cells at low Galactic 
latitude, |<5ga/| < 20°. In some cases we choose higher 
Galactic latitude cutoffs for comparison, as shown in Ta¬ 
bleland Figure© We also remove cells with high stel¬ 
lar density, or in very dusty regions, as described below. 
Although this does not explicitly avoid the Magellanic 
clouds, the high stellar density cutoffs remove their cen¬ 
tral regions. 

To examine the completeness of the galaxy sample, 
we have measured {V/Vmax) as a function of Ks,o rnag- 
nitude. The results are essentially constant between 
about 0.52 and 0.51 for 10.0 < A's.o < 13.7 and we em¬ 
ploy a conservative truncation of the galaxy sample at 
Ks,o = 13.5. 

To obtain the distribution function from counts-in- 
cells, we map the so urces onto a Hammer-Aitoff equal 
area projection fe.g.. IGalabretta fc Greisenll2002D : 


_ 180 ° / 2 

TT y 1 -I- cos(5cos(a/2) ’ 
x = 2'y cos^ sin(a/2), 

?/ = 7 sin 5. 

Here (a , 5) is the Galactic coordinate in radians and (x, 
y) is the projected coordinate. Figure©shows the sample 
of galaxy candidates with ATg.o < 13.5 on a square-root 

^ The Explanatory Supplement to the 
2MASS All Sky Data Release is available at 
http://www.ipac.caltech.edu/2mass/releases/allsky/doc for 

this and related information. 

^ The corrected magnitudes are denoted by a subscripted 0. 
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Fig. 1.— Grey-scale image of the density of 663,166 galaxy candidates with < 13.5 in cells of 1° x 1°. A square-root stretch having 
a contrast of 2 was applied with the darkest cells indicating the highest density. 


stretched grey scale with the darkest cells indicating the 
highest density. The extended dense region in the upper 
right hand (northwest) quadrant is the Shapley Super¬ 
cluster. This map contains 663,166 galaxy candidates in 
cells of 1° X l°on the sky. We then apply two filters. The 
first is a spatial filter that removes cells below a speci- 
hed Galactic latitude or on the boundary of the projec¬ 
tion. The second filter reduces the probability of confu¬ 
sion with stars by removing cells whose stellar density, 
Ust, is greater than 200 0 deg~^) and cells with 

objects that have A{Ka) > 0.05 l|Maller et al.ll200,^ . We 
experimented with a range of stellar density cutoffs and 
found nst < 10^'^ to be the optimum between leaving 
too few galaxies or including so many stars that a small 
percentage of them were close enough to each other that 
they could masquerade as galaxies. 

Having determined all the relevant cells, we simply 
count the number of galaxies in each to construct his¬ 
tograms of the distribution functions for various samples. 
These filtered samples include essentially the whole sky, 
four hemispherical sections, and four quadrant sections, 
along with a range of cell sizes, Galactic latitude cut¬ 
offs and magnitude limits. Values of N from the average 
projected density and by from the variance of counts in 
cells in equation 0 are used as a Hrst approximation to 
determine /v(iV) in equation Q) from the observations. 
Then we go a step further and Hnd values of N and b that 
minimize a fo observations. The differences in 
these values for a single histogram, and among similar 
histograms in different areas of the sky, give a measure of 
their global uncertainty. However values of cannot be 
used here for the usual probability estimates that a sam¬ 
ple represents a particular distribution function because 


the populations of nearby cells are often strongly corre¬ 
lated. Since these cells are not independent, the values 
of x^ are used just to optimize the htting. To facilitate 
this, we have combined bins w ith very small fv(N) i n 
the tails of the distribution (cf. iKenda.H V Stuarlill 

3. RESULTS 

To illustrate the relation between the predicted form of 
fv{N) in equation (^1 and the observed 2MASS galaxy 
distribution function, Figures[2lL-|2t plot least-squares fits 
of fv{N) to the observed histograms for square cells of 
hve different size cells on the sky. The areas under these 
histograms in the bottom panel of each hgure illustrate 

values of X = {fobs - ftheory)Nll^J, where Nyeu is 
the total number of cells evaluated. The 1° x 1° cells 
in FiguresEt andEli are representative; in FigureEJi, we 
have shifted the grid of cells by 0.5° on the sky. This shift 
slightly alters the detailed shape of the histogram for the 
probability that a cell has N galaxies, but leaves the best 
fit values of N and b almost unchanged. Thus the exact 
positioning of the grid is unimportant, which indicates 
that the galaxy distribution is statistically homogeneous. 

In the overall fits there are small differences between 
the observations and the theoretical distribution of equa¬ 
tion O- The largest differences, as measured by X; occur 
in the large N tails where there are relatively few cells. 
The peaks of the observed distribution function are usu¬ 
ally higher than those of equation O- Since the areas 
of the distributions have the same normalization, a slight 
excess for some values of N must be offset by a deficit for 
other values. This slight dehciency generally occurs be¬ 
tween the peak and high N tail of the distribution. There 
may be several contributions to these differences: the pe- 
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culiarity of the unique configuration of large clusters in 
the Universe, our emphasis on fitting the overall distribu¬ 
tion rather than giving greater weight to either the peak 
or tail, very small systematic effects in the 2MASS cat¬ 
alog and its analysis, or a real physical effect. If the last 
possibility is true, it might be caused by several effects 
including environmental interactions which tend to en¬ 
hance the luminosities of very overdense regions, mergers 
in relatively overdense regions which enhance the number 
of cells around the peak of the distribution, a remnant 
of initial conditions when galaxies started clustering, or 
segregation of more massive galaxies into smaller vol¬ 
umes. Since the overall differences are so small, it will 
be a considerable challenge to determine their causes. 

Figures -d show variations of these counts in cells. 
The upper two comparisons are for 1° x 1° cells (which we 
will now use for succeeding representative illustrations) 
over the entire sky, as before, but with magnitude cutoffs 
of Ksfi < 12.5 and Kgfi < 13.0. Examining the values 
of b for different magnitude cutoffs can describe whether 
magnitudes are partially correlated with position, as dis¬ 
cussed in Section ^ Comparison with Figure |21 shows 
that the excellence of the fit is not much affected by the 
magnitude cutoff, as we would expect for a statistically 
complete sample. Figures 0: and I3i show samples with 
the standard magnitude cutoffs but confined to Galactic 
latitudes more than 45° and 70° from the Galactic plane. 
Comparison with the 1° x 1° cells of Figures and 131, 
shows that the fits are as good and the values of h are 
the same within 2% for the 45° Galactic latitude cutoff 
and within 4% for the 70° cutoff. The 70° cutoff has 
larger fluctuations because it has only 20% of the num¬ 
ber of cells as the 45° cutoff and less than 10% of the 
cells in the 20° cutoff. Evidently galactic obscuration is 
not a significant problem above ±20° latitude for these 
distribution functions. 

The 2MASS sample is the first that is large enough 
for a precise determination of variations among different 
parts of the sky. Figures^^ - d give our results for divid¬ 
ing the sky into four independent quadrants. Each quad¬ 
rant contains between about 5,500 and 6,200 usable cells. 
Fluctuations around equation m are very small, though 
naturally they differ in detail. Among these quadrants 
the variance of N is 1.2% and of b is just 1.9%. For these 
quadrants divided into 4° x 4° cells, for which each quad¬ 
rant has between about 250 and 330 cells, the variance 
of N is 1.5% and of b is 2.6%. The smaller fluctuations 
among the larger cells, which are more representative of 
the total galaxy distribution, partly compensate for the 
smaller total number of such cells. This suggests that a 
very conservative estimate of the uncertainty in b is 3%. 
By averaging the two quadrants with the lowest values 
of b to give 0.485 and comparing with the averages of the 
two quadrants giving the largest value of 6 = 0.494, we 
would obtain an estimate of 2% for the uncertainty of 
b in cells of one square degree. For comparison, earlier 
estimates of b from the catalogs mentioned in Section ^ 
gave uncertainties of 5-10% for b. 

It is interesting to see the effects of a large cluster, in 
particular the Shapley Supercluster, on these statistics. 
On small scales it is obviously noticeable, as pointed out 
in the northwest quadrant of Figure ^ If it were impor¬ 
tant when averaged over the whole sky, or even over a 
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Fig. 3. — As Figure[3 except (a) and (b) indicate different mag¬ 
nitude cuts while (c) and (d) indicate different Galactic latitude 
cuts, all with 1° cells on the sky. The dotted lines in (c) and (d) 
indicate the best-fit model from the full sample with values of b 
and N given in Table [Tl 
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Fig. 4.— As Figure except the figures indicate different 
quadrants all with 1° cells on the sky. The dotted lines indicate 
the best-fit model from the full sample with values of b and N given 
in Table [T] 


quadrant, we would consider it a sufficient anomaly to 
question whether it could have formed by gravitational 
many-body clustering alone, or whether it needed to arise 
from large perturbations in the early universe. 

Figures |3t, b show the effects of removing the Shap- 
ley Supercluster from the all sky sample and from the 
northwest quadrant. This is effected by removing all cells 
within 10 degrees of {agai,Sgai) = (312.457°, 30.754°). 
Comparison with the same all sky samples in Figures[2t, 
d shows that removing the SSC reduces N slightly, as ex¬ 
pected, but reduces the value of b by less than 1%, and 
changes the histograms by less than the effects of the 
bin shift. On the smaller scale of a quadrant, however, 
there is less dilution and the effects of the SSC are more 
noticeable. The value of b is now reduced by 3%, com¬ 
pared with Figure 21 d, although the histograms remain 
very similar. 

A small fraction of anomalies in the 2MASS Galaxy 
(extended source) Catalog may be produced by contami¬ 
nation from point sources so close to each other that they 
masquerade as an extended source. For this to happen, 
their colors and magnitudes would also have to coincide 
with our criteria described in SectionElfor a galaxy. Fig¬ 
ures Et andEt, d illustrate this effect. In Figure Et we 
employ the less stringent criterion that all cells contain 
less than 10^ ® point sources, compared to limits of 10^'^ 
and lO^o in Figures Et and EH. This less stringent cri¬ 
terion induces slight fluctuations in the histograms but 
does not change the values of A^ or 5 to three signifi¬ 
cant figures. The more stringent criterion of Figure EH 
reduces both the number of cells and galaxies consider¬ 
ably (see Table n below) and has a 3% smaller value of 
b. This, and examination of Table E] suggest that the 
limit of 103-3 is a reasonable compromise. 

In Figure El we directly compare a sample from the 
2MASS Point Source Catalog (PSC, mainly stars) to a 
sample from the 2MASS XSC (mainly galaxies). To mit¬ 
igate the effects of the Galactic disc, the comparison uses 
a Galactic latitude cutoff of 70°. In the PSC, we first se¬ 
lected objects with valid photometry in all three bands, 
where Ks q < 13.5. Since there are many more objects 
in the PSC, we then randomly selected objects with a 
probability that would give approximately the same N 
as seen in the same area as the XSC. We also used the 
same stellar density requirements. The two distributions 
(XSC - solid line, PSC - dotted line) are clearly differ¬ 
ent with the PSC showing much less clustering (lower 
b). We also found that reducing the N for the PSC by 
magnitude selection produced a similar value of & « 0.05. 

The dependence of b on the scale and shape of the 
cells incorporates important information about the two- 
galaxy correlation function, as indicated by equation ®. 
Through the form of equation o, this dependence of 
b also incorporates much information about higher or¬ 
der correlations. This scale dependence can be used 
to test the applicability of computer simulations. In¬ 
deed a wide range of cold dark matter models do not 
agree with the scale de pendence of b in the IRAS Cat¬ 
alog llSheth et al ] IT^ . The increased precision of the 
2MASS Catalog now makes b{9) an even more stringent 
test of computer simulations. The b{9) dependence is 
shown in FigureEl where we have indicated conservative 
uncertainties of ±3% for each value of b. Typical CDM 
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Fig. 5. — As Figure HI except (a) and (b) indicate removing 
the Shapley Supercluster from the entire sky and the northwestern 
quadrant, while (c) and (d) indicate changing the stellar density 
at which cells are excluded, all with 1° size cells on the sky. The 
dotted lines indicate the best-fit model from the full sample with 
values of b and N given in Table [Tl 



Fig. 6.— Least-squares fits (continuous lines) of fyiN) to the 
observed histograms for likely galaxies (2MASS XSC, solid lines) 
and likely stars (2MASS PSC, dotted lines). 


models, especially those whose initial perturbation spec¬ 
tra produce large dark matter haloes that contain many 
galaxies, have too many voids and too many filaments at 
the present time (Baertschiger & Sylos-Labini 2004, pri¬ 
vate communication). Even if they agree approximately 
with the functional form of equation m, their excessive 
voids and filamentary structures would give significantly 
greater values of b on small or large scales than the obser¬ 
vations in Figured show. In Figured we also compare 
b{9) with values based on the two-galaxy angul ar corre¬ 
lation function, W{s). iLahav fc Saslawl lll992ll derived 
b[9) for square cells of size uj = 9 x 9 square degrees 

b{9) = 1-{1 + N sy^ (6) 

where W{s) = N is the expected number of 

galaxies projected into a cell, and is a coefficient 
evaluated numerically from J = ff 9^~'^ duji duj 2 = 
Cj 9^~^ as injTgtgjy^^Kihg^ l|1 96911 giving = 2.25 
for 7 = 1.8. fKOleretanH2flflT find 7 = 1.79 ±0.02 and 
So = 0.054±0.008 for the 2MASS catalog at < 13.5, 
from which we derive the dotted line in Figured This 
fit is less good at smaller 9 and has a = 33.6 for 5 
degrees of freedom. Fitting 7 and sq from 6(0), we find 
7 = 1.81 ± 0.04 and sq = 0.045 ± 0.003 (solid line in 
Figure ITll. with a = 3.0 , overlapping with the val¬ 
ues from iMa, Her et “^ ( 120 ?)^ . The distribution function, 
equation o, therefore provides a good method for mea¬ 
suring the angular correlation function. 

Tabledbelow gives a more detailed summary of results 
for a wider range of analyses than Figures dd The first 
section of Tabled examines samples with different stellar 
density limits in 1 x 1 ° cells over all the unobscured sky 
above ±20° Galactic latitude but with different magni¬ 
tude cutoffs. Values of N and b are determined in three 
different ways. In the first column of values, N is found 
exactly by averaging its value over all cells of each sam¬ 
ple, and 6 is found from the variance of counts in these 
cells using equation 0 - In the second column, N is 
found as before but 6 is found using a least-squares fit 
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Cell Length 9 (deg) 

Fig. 7.— Plot of b{9). The error bars represent conservative uncertain ties of di3% for each value of b. The predicted b{6) fequation IHl 
for the two-galaxy angular correlation function as measured for 2MASS in IMaller et all I2Q031 (dotted-line) and in this paper (solid-line). 


to equation O- In the third column, N and b are both 
obtained from a two-parameter least-squares fit to equa¬ 
tion O- In all cases, the fits are visually as good as those 
illustrated in Figures 1213 We have also given values of 
per degree of freedom, although as mentioned in Sec- 
tionOthese values cannot be used to calculate absolute or 
relative probabilities for comparisons among themselves 
or with other distributions. Moreover, these values 
tend to be dominated by the large N tail of the distri¬ 
butions. The values of b do not depend significantly on 
whether or not V is a fitting parameter or is directly 
determined from the counts in cells. The values of b de¬ 
termined from the variance using equation 0 tend to be 
high for small scales but agree well with the values from 
the fits for larger cells which are more representatiave 
(even though there are fewer of them). The remaining 
sections of Tablenexamine the effects of varying Galactic 
latitude cutoffs, varying the size of the cell, and varying 
the regions of the sky to which the GQED is fit. 

In Table EJ we compare the effects of sampling the 
galaxy catalog by magnitude cuts with sampling by a 
random selection that yields approximately the same 
number of galaxies at each magnitude cut. For these 
comparisons, we have used the usual stellar density cut¬ 
off {ust < 10^-^), All-sky region, and Galactic latitude 
cutoff (above ±20°). The value of 6 in a subsample se¬ 
lected randomly with a probabilit y p is related to the 
value of the parent sample by (e.g., ISa,s1a,wll2?innll 


(1 bsubsaraple) — 


(1 bparen9) 


^ P){‘^ ^paren9bparent 


( 7 ) 


The values of b in Table El for the random sampling of 
the 13.5 magnitude limited galaxy catalog compared with 
the values of b in Table E] confirm equation (TJ . The 
pairs of rows in Table El give observed parameters for a 
magnitude cut and for a corresponding Poisson selected 
subsample. This analysis is for cells of 1 x 1°, 2 x 2°, 
and 4x4°. Figure El shows the differences between b 
measured for the magnitude cut and b for an equivalent 
randomly selected subsample. 

The magnitude selected subsample always has a higher 
observed value of b than the corresponding randomly se¬ 
lected subsample. This effect appears to be mainly ge¬ 
ometric. As one goes to brighter subsamples, their av¬ 
erage volume decreases provided the distribution is sta¬ 
tistically homogeneous. This affects W{s) by increasing 
Sq. Assuming the average radial distance sampled follows 
D{m) oc 10™/^, and with the usual scal ing argument, 
W{s) (X W{s D{m))/D{m) llPeebleiJII993|l . 

so.m, (8) 

Combining equations El and El we can predict b{9, m) 
based on W{s) at Kgfi = 13.5. In Figure El we show the 
measured b minus this predi cted b. Using the a ngular 
correlation, IF(s, 13.5), from iMaller et all l)2005ll . there 
are some small discrepancies for small 9. With our best- 
fit W{s, 13.5), there are no discrepancies. 

4. SUMMARY AND DISCUSSION 

Our main result is that at a level of precision 2-3 times 
better than possible with previous catalogs, the detailed 
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Fig. 8.— Plot of the differences between values of b measured from magnitude limited subsamples and from randomly selected subsamples 
of equal size as a function of magnitude for three cell sizes. Small shifts in plotted magnitude are made so that individual data-points and 
error bars can be distinguished. The error bars represent conservative uncertainties scaled by N to ±3% of b for Ks^q = 13.5. 


form of the galaxy distribution in the 2MASS Catalog 
agrees with its prediction from gravitational cosmologi¬ 
cal many-body theory. The value of the clustering pa¬ 
rameter b in equation Q appears to be 6 = 0.867±0.026 
on large scales in the present universe with z < 0.1. This 
is in accordance with the value of 0.83 ± 0.05 for the 
much smaller but also h ighly complete three-dimensiona l 
Pisces-Perseus sample llSaslaw fc Haaue-CoDilahlll998jl . 
Tabled shows that as the cell size increases the values of 
b deduced from equation 0 agree better with the values 
in the fits. This is essentially because larger cells are 
more representative of the global clustering, and their 
variance is more accurately determined. These values of 
b also agree with N-body simulatio ns at the pres ent time 
in rio = 1 models (summarized in ISa,s1a.wll2(lflnll . Agree¬ 
ment with gravitational many-body theory need not have 
occurred at this level, especially since there are at least 
four major processes that could have destroyed it. The 
reasons that the agreement survives may be quite in¬ 
structive. 

First, a non-uniform distribution of dark gravitating 
matter that differs from the luminous galaxy distribu¬ 
tion could significantly modify the observed galaxy dis¬ 
tribution, in disagreement with the form of equation O- 
If dark haloes contain many galaxies, the resulting voids 
and filaments could exceed those observed. The value of b 
in such models would be too high. Moreover models with 
overlapping individual dark matter haloes would soften 
the gravitational potential significantly and their result¬ 


ing peculiar velocity dist ribution function would also dis¬ 
agree with observations l|T;enng fc Sa,s1a,wll2f)n4jl . All this 
indicates that when galaxies cluster most of the dark 
matter is located in the haloes of individual galaxies. 

Second, if the initial perturbations that stimulated 
galaxy clustering have survived into the present non¬ 
linear regime, they could produce departures from the 
form of equation O- This would be especially true if 
the initial conditions were strongly correlated or anti¬ 
correlated over large scales and have not had time to 
relax through a sequence of quasi-equilibrium states 
llSaslawl 1200(111 . For example, an initial power law per¬ 
turbation spectrum would have to have an expone nt be¬ 
tween — 1 < n < 1 for such relaxation l|Itohlll9^ 1. Of 
course the initial power spectrum may not have been 
a simple power law; then this question becomes much 
more complicated and has not yet been systematically 
explored. 

Third, there is considerable evidence, especially from 
the Hubble Deep Fields, that merging has played a sig¬ 
nificant role in galaxy formation which extends into the 
present. Merging tends to modify the form of fv(N) 
slight ly and decrease the value of b llFang fc Saslawl 
[1993). The effects of merger induced starbursts modify 
the luminosity function, and complicate the selection of 
galaxies in the distribution function, probably depending 
on the local environment. 

Evidently mergers do not destroy the distribution func¬ 
tions in Section|21 There may be several reasons for this. 
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Fig. 9. — Plot of the differences between values of b measured from magnitude limited subsamples and predictions for b from equations 
|H|_andJSJasa function of magnitude for three cell sizes. We use the two-galaxy angular correlation function as measured for 2MASS in 
|MaJler_et_aJJ J2003) (dotted error-bars) and in this paper (solid error-bars). Small shifts in plotted magnitude are made so that individual 
data-points and error bars can be distinguished. The error bars represent conservative uncertainties scaled by N to ±3% of b for = 13.5. 


Although individual mergers may be dramatic, only a 
small percentage of galaxies merge at any given time so 
they do not have much effect on the overall statistics. 
When galaxies merge, they generally do so near the cen¬ 
ter of mass with their nearest neighbor. These centers 
of mass may have a distribution which approximately 
follows that of the galaxies themselves. Moreover in¬ 
tervals between mergers are often long compared with 
the local orbital relaxation time, so the merged galaxy 
has time to acquire the quasi-equilibrium distribution of 
equation O- 

Fourth, although the distribution of equation o is for 
systems in which all galaxies have t he same mass, N-body 
simulations (e.g.. lItoh et al.llT99,‘lD show that a range of 
masses generally enhances clustering and increases the 
value of b compared with single mass cases. The effect 
is relatively small, typically ~10%, because much of the 
later clustering involves the collective interaction of indi¬ 
vidual galaxies with groups of galaxies, or of groups with 
other groups, for which the masses of individual galaxies 
are less relevant. For a wide range of masses, there is 
a tendency for massive galaxies to acquire less massive 
galaxies as satellites. How much this affects the distri¬ 
bution function depends on the relative total numbers of 
low and high mass galaxies and, for the observations, on 
detailed models of mass-luminosity galaxy selection. 

While it is conceivable that these four, and perhaps 
other, effects might conspire to leave the distribution 


function of equation o invariant, it seems more likely 
to us that each effect is individually small. Nonetheless 
they may be present to some degree since the fits to equa¬ 
tion o are not perfect. In particular, there seem to be 
slight systematic tendencies for the observed peaks to ex¬ 
ceed the theoretical peaks, for a small observed excess in 
the large N tail, and for a mild observed deficit between 
the peak and this tail. This holds for many samples, but 
not all. It seems somewhat more prominent in samples 
with greater stellar contamination, but might also be pro¬ 
duced by some of the effects mentioned previously. At 
present, it seems that tracking down these small discrep¬ 
ancies will require a combination of even more extensive 
observations and more detailed models. 

We have also found (Table |3) that values of b for sub¬ 
samples limited by brighter apparent magnitude cutoffs 
are higher than b would be for a subsample selected ran¬ 
domly with the sa me probability. This wa s not evident in 
an earlier analysis llSa,s1a,w fc Cra,ne|ll 99ll) of the Zwicky 
catalog, probably because this earlier catalog contains 
only about 7% as many galaxies as 2MASS and differ¬ 
ences of < I mag were explored. This effect can be ex¬ 
plained almost entirely by the smaller volume of sub¬ 
samples with brighter magnitude cutoffs, as Figure El 
shows. From equation ®, the spatially integrated effect 
of luminosity clustering as observed in the luminosity 
dependence of the two-point correlation function (e.g. . 
iBenoist et alJII99^ iMo et al1lI994 iNorberg et alJl2?)?)ll) 






































is therefore small. Similarly, the good fits of our magni¬ 
tude limited subsamples to equation o, which includes 
all the higher order correlation functions and moments of 
the distribution, show that luminosity clustering is small 
when averaged over large statistically homogeneous vol¬ 
umes containing many clusters. Physically, this would 
be expected if most of the clustering we re produced by 
collective gravitational interactions (e.g., lSaslawll200nll . 

With further information about other properties of the 
2MASS galaxies, distribution functions will make possi¬ 
ble new insights into relative clusterings based on mor¬ 
phological types, colors, presence of satellite and tidally 
interacting galaxies, etc. This will contribute to a wealth 
of constraints on galaxy formation and evolution in more 
detailed models. 


This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Pro¬ 
cessing and Analysis Center/California Institute of Tech¬ 
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Administration and the National Science Foundation. 
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advice about the 2MASS survey and to the anonymous 
referee for especially helpful comments. GRS acknowl¬ 
edges the receipt of an Achievement Reward for College 
Scientists fellowship. Partial support was also provided 
by the F.H. Levinson Fund of the Peninsula Community 
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Table 1. GQED Fits to 2MASS Data 


Max 

Max 


Min 


Total 

Total 

From Variance 


Eitting b 

Fitting A & 6 

Kso 

log{nst) 

Region 

\^Gal\ 

8 

Galaxies 

Cells 

N 

b 

xVDof 

N 

b 

xVDof 

N 

b 

xVDof 

For Varying rist, KsO Cutoffs 

11.00 3.0 All 20 

1.00 

6717 

15199 

0.44 

0.228 

26.4/ 7 

0.44 

0.203 

14.9/ 6 

0.44 

0.200 

14.2/ 5 

11.00 

3.3 

All 

20 

1.00 

10940 

23348 

0.47 

0.264 

84.5/ 8 

0.47 

0.223 

37.6/ 7 

0.46 

0.218 

34.8/ 6 

11.00 

3.6 

All 

20 

1.00 

11421 

24171 

0.47 

0.274 

116.0/ 8 

0.47 

0.223 

42.8/ 7 

0.46 

0.217 

39.0/ 6 

11.50 

3.0 

All 

20 

1.00 

13655 

15199 

0.90 

0.283 

55.9/ 9 

0.90 

0.248 

23.9/ 8 

0.89 

0.245 

22.3/ 7 

11.50 

3.3 

All 

20 

1.00 

22252 

23348 

0.95 

0.329 

250.8/12 

0.95 

0.266 

87.3/11 

0.92 

0.257 

76.7/10 

11.50 

3.6 

All 

20 

1.00 

23131 

24171 

0.96 

0.333 

288.8/12 

0.96 

0.265 

91.7/11 

0.93 

0.256 

79.3/10 

12.00 

3.0 

All 

20 

1.00 

27872 

15199 

1.83 

0.340 

158.5/13 

1.83 

0.291 

72.9/12 

1.79 

0.284 

66.6/11 

12.00 

3.3 

All 

20 

1.00 

45385 

23348 

1.94 

0.388 

479.2/17 

1.94 

0.317 

166.8/16 

1.88 

0.305 

144.5/15 

12.00 

3.6 

All 

20 

1.00 

47145 

24171 

1.95 

0.389 

505.4/17 

1.95 

0.317 

169.6/16 

1.88 

0.305 

145.6/15 

12.50 

3.0 

All 

20 

1.00 

58035 

15199 

3.82 

0.398 

193.0/20 

3.82 

0.355 

94.8/19 

3.75 

0.349 

88.0/18 

12.50 

3.3 

All 

20 

1.00 

94467 

23348 

4.05 

0.453 

731.8/24 

4.05 

0.382 

243.3/23 

3.92 

0.368 

209.6/22 

12.50 

3.6 

All 

20 

1.00 

97926 

24171 

4.05 

0.453 

751.8/24 

4.05 

0.382 

247.3/23 

3.92 

0.369 

212.5/22 

12.75 

3.0 

All 

20 

1.00 

84186 

15199 

5.54 

0.432 

204.7/25 

5.54 

0.392 

98.8/24 

5.46 

0.386 

92.7/23 

12.75 

3.3 

All 

20 

1.00 

136610 

23348 

5.85 

0.485 

836.6/29 

5.85 

0.417 

264.1/28 

5.69 

0.405 

232.3/27 

12.75 

3.6 

All 

20 

1.00 

141599 

24171 

5.86 

0.485 

851.7/29 

5.86 

0.417 

268.0/28 

5.70 

0.406 

235.6/27 

13.00 

3.0 

All 

20 

1.00 

122608 

15199 

8.07 

0.461 

271.4/31 

8.07 

0.422 

143.9/30 

7.95 

0.416 

136.1/29 

13.00 

3.3 

All 

20 

1.00 

197738 

23348 

8.47 

0.514 

1108.0/37 

8.47 

0.446 

364.2/36 

8.24 

0.433 

326.3/35 

13.00 

3.6 

All 

20 

1.00 

204826 

24171 

8.47 

0.514 

1111.1/37 

8.47 

0.447 

366.9/36 

8.25 

0.434 

328.1/35 

13.25 

3.0 

All 

20 

1.00 

178009 

15199 

11.71 

0.492 

264.8/38 

11.71 

0.456 

138.3/37 

11.55 

0.450 

129.6/36 

13.25 

3.3 

All 

20 

1.00 

285786 

23348 

12.24 

0.540 

1073.2/45 

12.24 

0.479 

359.9/44 

11.94 

0.467 

320.5/43 

13.25 

3.6 

All 

20 

1.00 

296039 

24171 

12.25 

0.540 

1082.8/46 

12.25 

0.479 

369.2/45 

11.95 

0.467 

328.2/44 

13.50 

3.0 

All 

20 

1.00 

258815 

15199 

17.03 

0.520 

308.5/47 

17.03 

0.485 

168.0/46 

16.83 

0.480 

159.9/45 

13.50 

3.1 

All 

20 

1.00 

333937 

19170 

17.42 

0.541 

617.0/52 

17.42 

0.495 

275.6/51 

17.14 

0.487 

256.9/50 

13.50 

3.2 

All 

20 

1.00 

382096 

21652 

17.65 

0.554 

831.0/54 

17.65 

0.504 

346.7/53 

17.32 

0.495 

318.9/52 

IS.SO'^ 

3.3 

All 

20 

1.00 

413339 

23348 

17.70 

0.565 

1126.1/54 

17.70 

0.507 

397.7/53 

17.34 

0.496 

359.9/52 

13.50 

3.4 

All 

20 

1.00 

422901 

23894 

17.70 

0.564 

1130.7/55 

17.70 

0.507 

404.5/54 

17.33 

0.497 

365.9/53 

13.50 

3.5 

All 

20 

1.00 

427446 

24136 

17.71 

0.563 

1142.7/56 

17.71 

0.507 

408.4/55 

17.34 

0.496 

369.0/54 

13.50 

3.6 

All 

20 

1.00 

428109 

24171 

17.71 

0.563 

1137.5/56 

17.71 

0.507 

406.5/55 

17.34 

0.496 

367.2/54 

13.50 

4.0 

All 

20 

1.00 

428578 

24194 

17.71 

0.563 

1142.0/56 

17.71 

0.507 

408.1/55 

17.34 

0.497 

368.5/54 

For Varying Galactic Latitude Cutoffs 
13.50^ 3.3 All 20 1.00 

413339 

23348 

17.70 

0.565 

1126.1/54 

17.70 

0.507 

397.7/53 

17.34 

0.496 

359.9/52 

13.50 

3.3 

All 

45 

1.00 

203044 

11439 

17.75 

0.550 

452.0/48 

17.75 

0.499 

181.3/47 

17.42 

0.490 

165.2/46 

13.50 

3.3 

All 

70 

1.00 

39548 

2223 

17.79 

0.555 

187.1/32 

17.79 

0.488 

92.6/31 

17.29 

0.472 

85.7/30 

For Varying 9 
13.50 3.3 

All 

20 

0.25 

439754 

412101 

1.07 

0.258 

2266.4/19 

1.07 

0.217 

975.6/18 

1.05 

0.212 

897.9/17 

13.50 

3.3 

All 

20 

0.50 

430439 

98048 

4.39 

0.406 

2069.1/31 

4.39 

0.349 

774.9/30 

4.29 

0.338 

688.3/29 

IS.SO'^ 

3.3 

All 

20 

1.00 

413339 

23348 

17.70 

0.565 

1126.1/54 

17.70 

0.507 

397.7/53 

17.34 

0.496 

359.9/52 

IS.SO'* 

3.3 

All 

20 

2.00 

384240 

5409 

71.04 

0.697 

272.8/93 

71.04 

0.662 

138.0/92 

69.81 

0.655 

126.1/91 

13.50^ 

3.3 

All 

20 

4.00 

330298 

1159 

284.99 

0.793 

91.7/43 

284.99 

0.771 

61.7/42 

283.80 

0.770 

61.4/41 

13.50 

3.3 

All 

20 

6.00 

280916 

435 

645.78 

0.835 

13.4/17 

645.78 

0.828 

11.7/16 

643.56 

0.828 

11.6/15 

13.50 

3.3 

All 

20 

8.00 

231854 

201 

1153.50 

0.858 

3.5/ 7 

1153.50 

0.867 

2.6/ 6 

1154.00 

0.867 

2.6/ 5 

For Varying Regions of the 
13.50^ 3.3 All 

Sky Using 8 = i 
20 1.00 

A Cells 
413339 

23348 

17.70 

0.565 

1126.1/54 

17.70 

0.507 

397.7/53 

17.34 

0.496 

359.9/52 

13.50 

3.3 

All 

20 

l.OO*’ 

413308 

23365 

17.69 

0.562 

995.5/55 

17.69 

0.508 

352.6/54 

17.33 

0.498 

314.7/53 

13.50 

3.3 

AIF 

20 

1.00 

404863 

23050 

17.56 

0.554 

951.0/54 

17.56 

0.502 

374.4/53 

17.22 

0.492 

341.0/52 

13.50 

3.3 

N 

20 

1.00 

217041 

12220 

17.76 

0.573 

681.0/49 

17.76 

0.512 

221.7/48 

17.39 

0.502 

202.0/47 

13.50 

3.3 

S 

20 

1.00 

196298 

11128 

17.64 

0.553 

501.7/48 

17.64 

0.500 

231.9/47 

17.24 

0.489 

210.5/46 

13.50 

3.3 

E 

20 

1.00 

201571 

11465 

17.58 

0.561 

523.1/48 

17.58 

0.509 

221.0/47 

17.21 

0.498 

202.8/46 

13.50 

3.3 

W 

20 

1.00 

210012 

11786 

17.82 

0.566 

643.1/49 

17.82 

0.504 

203.8/48 

17.43 

0.493 

181.6/47 

13.50 

3.3 

NW 

20 

1.00 

111595 

6219 

17.97 

0.562 

327.6/42 

17.97 

0.501 

103.4/41 

17.65 

0.494 

97.1/40 

13.50 

3.3 

NW‘= 

20 

1.00 

103119 

5921 

17.42 

0.520 

327.6/39 

17.42 

0.483 

76.6/38 

17.26 

0.479 

47.7/37 

13.50 

3.3 

NE 

20 

1.00 

104720 

5958 

17.58 

0.584 

370.9/41 

17.58 

0.523 

138.4/40 

17.11 

0.509 

124.3/39 






Table 1. GQED Fits to 2MASS Data— Continued 


Max 

Max 


Min 


Total 

Total 

From Variance 


Fitting b 



Fitting N b 

Kso 

log{nst) 

Region 

\Saal\ 

9 

Galaxies 

Cells 

N 

b 

X^/Dof 

N 

b 

X^/Dof 

N 

b 

xVDof 

13.50 

3.3 

SE 

20 

1.00 

96851 

5507 

17.59 

0.532 

198.5/40 

17.59 

0.493 

131.1/39 

17.28 

0.484 

124.4/38 

13.50 3.3 SW 20 1.00 

For Varying Regions of the Sky Using 0 = 4° 

98417 

Cells 

5567 

17.68 

0.570 

362.7/42 

17.68 

0.505 

142.0/41 

17.16 

0.490 

124.0/40 

13.50^ 

3.3 

All 

20 

4.00 

330298 

1159 

284.99 

0.793 

91.7/43 

284.99 

0.771 

61.7/42 

283.81 

0.770 

61.4/41 

13.50 

3.3 

NW 

20 

4.00 

95024 

328 

289.71 

0.797 

1Q.2I12 

289.71 

0.776 

11.3/11 

286.57 

0.773 

10.8/10 

13.50 

3.3 

NE 

20 

4.00 

83880 

294 

285.31 

0.817 

8.1/11 

285.31 

0.802 

5.4/10 

282.30 

0.799 

5.0/ 9 

13.50 

3.3 

SE 

20 

4.00 

71567 

253 

282.87 

0.756 

2.7/ 9 

282.87 

0.747 

2.4/ 8 

279.20 

0.742 

1.7/ 7 

13.50 

3.3 

SW 

20 

4.00 

73875 

264 

279.83 

0.783 

14.8/10 

279.83 

0.771 

13.7/ 9 

275.24 

0.765 

12.7/ 8 


Note. — The units of Ust are deg The units of SQai and 0 are deg. 


This row is a parent sample for randomly selected subsamples in Table 
The Galactic longitude bins were shifted by 0.5° for this fit. 

Gells containing galaxies within 10 degrees the Shapley Supercluster were excluded for this fit. 









Table 2. Comparison of Magnitude Selected Subsamples with Ran¬ 
domly Selected Subsamples 


Max 

Kso 

Poisson 

Prob. 

e 

Total 

Galaxies 

Total 

From Variance 


Fitting b 

Fitting N k b 

Cells 

N 

b 

xVDof 

N 

b 

xVDof 

N 

b 

xVBoi 

For Varying 

1 Cutoffs Compared to Random Selection from KsO = 13.50 Cutoff Using 6 = V 

" Cells 




11.00 

1.0000 

1.00 

10940 

23348 

0.47 

0.264 

84.5/ 8 

0.47 

0.223 

37.6/ 7 

0.46 

0.218 

34.8/ 6 

13.50 

0.0265 

1.00 

10917 

23348 

0.47 

0.049 

12.1/ 5 

0.47 

0.044 

11.2/ 4 

0.47 

0.043 

11.1/ 3 

11.50 

1.0000 

1.00 

22252 

23348 

0.95 

0.329 

250.8/12 

0.95 

0.266 

87.3/11 

0.92 

0.257 

76.7/10 

13.50 

0.0538 

1.00 

22343 

23348 

0.96 

0.091 

13.3/ 8 

0.96 

0.085 

11.4/ 7 

0.95 

0.084 

11.3/ 6 

12.00 

1.0000 

1.00 

45385 

23348 

1.94 

0.388 

479.2/17 

1.94 

0.317 

166.8/16 

1.88 

0.305 

144.5/15 

13.50 

0.1100 

1.00 

45530 

23348 

1.95 

0.188 

87.5/12 

1.95 

0.158 

45.2/11 

1.93 

0.157 

43.2/10 

12.50 

1.0000 

1.00 

94467 

23348 

4.05 

0.453 

731.8/24 

4.05 

0.382 

243.3/23 

3.92 

0.368 

209.6/22 

13.50 

0.2290 

1.00 

94631 

23348 

4.05 

0.292 

278.1/19 

4.05 

0.246 

142.6/18 

4.00 

0.241 

132.2/17 

13.00 

1.0000 

1.00 

197738 

23348 

8.47 

0.514 

1108.0/37 

8.47 

0.446 

364.2/36 

8.24 

0.433 

326.3/35 

13.50 

0.4780 

1.00 

197645 

23348 

8.47 

0.425 

604.3/33 

8.47 

0.371 

262.7/32 

8.31 

0.362 

239.9/31 

For Varying 

1 Cutoffs Compared to Random Selection from KsO = 13.50 Cutoff Using 0 = 2' 

" Cells 




11.00 

1.0000 

2.00 

10111 

5409 

1.87 

0.390 

49.3/12 

1.87 

0.352 

29.6/11 

1.82 

0.344 

27.3/10 

13.50 

0.0263 

2.00 

10279 

5409 

1.90 

0.112 

10.3/ 8 

1.90 

0.102 

9.2/ 7 

1.90 

0.102 

9.1/ 6 

11.50 

1.0000 

2.00 

20585 

5409 

3.81 

0.465 

85.9/18 

3.81 

0.418 

38.2/17 

3.71 

0.409 

34.5/16 

13.50 

0.0536 

2.00 

20859 

5409 

3.86 

0.179 

25.3/13 

3.86 

0.154 

16.9/12 

3.84 

0.152 

16.5/11 

12.00 

1.0000 

2.00 

42092 

5409 

7.78 

0.537 

205.3/27 

7.78 

0.481 

95.8/26 

7.53 

0.467 

85.7/25 

13.50 

0.1100 

2.00 

42268 

5409 

7.81 

0.309 

34.7/20 

7.81 

0.281 

20.9/19 

7.76 

0.279 

20.0/18 

12.50 

1.0000 

2.00 

87685 

5409 

16.21 

0.603 

297.3/41 

16.21 

0.547 

127.1/40 

15.69 

0.532 

111.0/39 

13.50 

0.2280 

2.00 

87784 

5409 

16.23 

0.450 

162.2/34 

15.23 

0.401 

84.5/33 

16.02 

0.394 

79.7/32 

13.00 

1.0000 

2.00 

183757 

5409 

33.97 

0.656 

342.5/61 

33.97 

0.610 

164.9/60 

33.15 

0.599 

150.2/59 

13.50 

0.4780 

2.00 

183759 

5409 

33.97 

0.579 

203.4/56 

33.97 

0.539 

114.1/55 

33.43 

0.532 

104.5/54 

For Varying 

1 Cutoffs Compared to Random Selection from KsO = 13.50 Cutoff Using 0 = 4' 

" Cells 




11.00 

1.0000 

4.00 

8509 

1159 

7.34 

0.546 

37.3/17 

7.34 

0.504 

22.7/16 

7.19 

0.496 

21.9/15 

13.50 

0.0258 

4.00 

8517 

1159 

7.35 

0.217 

12.4/14 

7.35 

0.203 

11.8/13 

7.32 

0.202 

11.7/12 

11.50 

1.0000 

4.00 

17430 

1159 

15.04 

0.610 

55.1/24 

15.04 

0.575 

40.3/23 

14.63 

0.564 

38.2/22 

13.50 

0.0528 

4.00 

17385 

1159 

15.00 

0.325 

23.2/20 

15.00 

0.304 

21.2/19 

14.93 

0.301 

21.0/18 

12.00 

1.0000 

4.00 

35665 

1159 

30.77 

0.674 

96.1/35 

30.77 

0.633 

64.0/34 

29.95 

0.622 

60.9/33 

13.50 

0.1080 

4.00 

35619 

1159 

30.73 

0.449 

53.2/29 

30.73 

0.424 

49.1/28 

30.47 

0.420 

48.3/27 

12.50 

1.0000 

4.00 

74812 

1159 

64.55 

0.725 

78.9/35 

64.55 

0.689 

40.9/34 

63.61 

0.683 

39.6/33 

13.50 

0.2260 

4.00 

74913 

1159 

64.64 

0.591 

57.8/33 

64.64 

0.557 

41.7/32 

64.45 

0.555 

41.6/31 

13.00 

1.0000 

4.00 

157519 

1159 

135.91 

0.764 

77.1/39 

135.91 

0.735 

40.4/38 

134.46 

0.731 

39.4/37 

13.50 

0.4770 

4.00 

157197 

1159 

135.63 

0.708 

66.5/39 

135.63 

0.680 

43.8/38 

134.58 

0.677 

42.9/37 


Note, 


All fits were to All-sky data where Ust < 3.3 deg ^ and I^Cail > 20 deg. The units of d are deg. 



